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Abstract. We derive the energy levels associated with the even-parity wave functions 
of the harmonic oscillator with an additional delta-function potential at the origin. 
Our results bring to the attention of students a non-trivial and analytical example 
of a modification of the usual harmonic oscillator potential, with emphasis on the 
modification of the boundary conditions at the origin. This problem calls the attention 
of the students to an inaccurate statement in quantum mechanics textbooks often found 
in the context of solution of the harmonic oscillator problem. 
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1. Introduction 

Every single book on quantum mechanics gives the solution of the harmonic oscillator 
problem. The reasons for that are, at least, two: (i) it is a simple problem, amenable 
to different methods of solution, such as the Frobenius method for solving differential 
equations [U [2], and the algebraic method leading to the introduction of creation and 
annihilation operators [3]. This problem has therefore a natural pedagogical value; (ii) 
the system itself has immense applications in different fields of physics and chemistry 
[U |5] and it will appear time and time again in the scientific life of a physicist. 

Another problem often found in quantum mechanics textbooks is the calculation 
of the bound state (negative energy) of the potential V(x) = a8(x), with a < 0. The 
latter example is instructive for the student because it cracks down the miss-conception 
that the continuity of the wave function and its first derivative at an interface is the 
only possible boundary condition in quantum problems. In fact, for this potential, the 
first derivative of the wave function is discontinuous at x — 0. To see the origin of this 
result, let us write the Schrodinger equation as 

Integrating Eq. ([I]) in an infinitesimal region around x = we get [3J 



lim 



e->0 dx 



^(0) = # v£(o) - ^<(0) = ^(0), (2) 
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where i/j^ (x) are the wave functions in each side of the delta potential and we have used 
Newton's notation for derivatives where the primes over functions denote the order of 
the derivative of that function. It is clear from Eq. ^ that the first derivative of the 
wave function is discontinuous. An explicitly calculation gives the eigenstate of the 
system in the form 

ijj{x) = v^e- KN , (3) 

which has a kink at the origin as well as a characteristic length scale given by k = ma/ti 2 . 
The eigenvalue associated with the wave function @ is E = -h 2 K 2 /(2m). When E > 0, 
the system has only scattering states for both positive and negative values of a. We 
now superimpose a harmonic potential potential on the already present delta-function 
potential. Due to the confining harmonic potential, the new system only has bound 
states, no mater the sign of a and E. Thus, the problem we want to address is the 
calculation of the eigenstates and eigenvalues of the potential 

777 

V(x) = -u¥ + a6(x) . (4) 

As we will see, this problem has both "trivial" and non-trivial solutions. 
Furthermore, it allows a little excursion into the world of special functions. Indeed, 
special functions play a prominent role in theoretical physics, to a point that the famous 
Handbook of Mathematical Functions, by Milton Abramowitz and Irene Stegun [6], would 
be one of the three texts (together with the Bible and Shakespeare complete works) 
Michal Berry would take with him to a desert island [7j. In a time where symbolic 
computational software is becoming more and more the source of mathematical data, 
we hope with this problem to show that everything we need can be found the good old 
text of Milton Abramowitz and Irene Stegun [6] . 

In addition, this problem will call the attention of students to an inaccurate 
statement in quantum mechanics textbooks often found in the context of the solution 
of the harmonic oscillator problem. Our approach is pedagogical, in the sense that 
illuminates the role of boundary conditions imposed on the even-parity wave functions 
by the 5— function potential. (We note that after the submission of this article, the work 
by Busch et al. [8] was brought to our attention; see note added at the end of article.) 



2. Solution of the harmonic oscillator with a delta-function 

The Hamiltonian of the system is 

H 2 d 2 m 

H = ~^ + T + "*(*>• («) 
Following tradition, we introduce dimensionless variables using the intrinsic length scale 
of the problem = h/{muj). Making the substitution y = x/ao in Eq. ^ we can write 
the Schrodinger equation as 

+ (2e _ jfMy) + 2g5{y)^{y) = , (6) 
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where e = malE/h 2 and g = aaom/h 2 . With g — 0, Eq. (|6| is recognized as the 
Weber- Hermite differential equation p]. In quantum mechanics textbooks, the solution 
of Eq. ^ with g = proceeds by making the substitution 

iKv) = e-y 2 ' 2 w{y) ■ (7) 

At same time, it is a common practice to write 2e = 2v + 1, where v is a real number. 
This allow us to transform Eq. ^ into 

w" — 2yw' + 2uw — 2g5(y)w = , (8) 

which is Hermite's differential equation when g = 0; a further substitution, 2 = y 2 , 
transforms the Eq.([8]) into the Kummer's equation, 

zw" + {b — z)w' — auw = with a — — | and b = | , (9) 

which, obviously, has two linearly independent solutions: the confluent hypergeometric 
functions M(a,b,z) and U(a,b,z); these functions also known as Kummer's functions 
(the latter solution is sometimes referred as Tricomi's function) [fj Thus, the general 
solution of Eq. ^ is 

w(z)=A v M(-%,l,z)+B v U(-%,±,z), (10) 

where A u and B u are arbitrary complex constants and v is an arbitrary real number. 
The U{— | , 2> 2) function can also be written in terms the functions M(a, b, z) as [9la] 

[ "" ' tr(l)r(i-l) V * r(l)r(-|) f 



It is important to note that Eq. (11) is not a linear combination of two M(a,b,z) 



functions. Using Eq. (11), it is possible to show that if v is either zero or a positive 



integer number (denoted by n), the solution in Eq. (10) can put in the following form 
(* > 0): 1 

J M(-|, \,z) for n even 

10 W 00 I ^M{\ - \,\,z) for n odd, ^ 

or more compactly [Sl6,c], it can be written as w(z) oc H n (y/z), with H n (z) the Hermite 
polynomial of order n; further more the product w(z)e~^ converges for all z. Then, the 
full wave function has the usual form 

with the corresponding eigenvalues being e = n + \. What we have detailed 
above condenses the typical solution of the quantum harmonic oscillator using special 
functions. 

X These functions may also be referred as the confluent hypergeometric function of the first and second 
kind, with the notation M(a, b, z) —1 Fi(a; b; z) and U(a, b, z) = ^ a F (a; 1 + a — b; — 1/z). All these 
notations can be found when using computational methods and software. 



§ This is easy to see from Eq. (11) recalling that the Gamma function diverges at negative integer 
values. 
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We now move to the solution of the quantum harmonic oscillator with a ^-function 
potential at the origin. To that end, we have to review few properties of the M(a, b, z) 
functions. For non-integer values of a and b, the M(a, b, z) function is a convergent 
series for all finite given z Pa], but diverges for z — > +00 as ^d] 

M(a, 6, z) = S^eV^[l + 0(\z\-% (13) 



In terms of the original variable y 2 , the function (13) diverges as e y2 , which implies that 
tp{y) also diverges at infinity as e y ^ 2 . Thus, ip(y) is not, in general, an acceptable wave 
function. 

We noted at the introduction to this article that it is many times referred 
(erroneously) in most standard textbooks on quantum mechanics that the only 
mathematical solutions of the harmonic oscillator differential equation that does not 
blow up when y — > +00 are those having v either zero or a positive integer. 

On the contrary however, the function U(— |, \ 1 y 2 ) with a non-integer v does not 
blow up as e y2 when y — >■ +00. Indeed, it is easy to see [using Eqs.(ll) and (13)] that 
U(—u/2, 1/2, y 2 ) -» y v [9,e] as y ->■ +00. Thus, the function 

^( y ) = Ae-h 2 u(-$,k,y 2 ), (14) 

could in principle be an acceptable wave function for any value of z/, since it is both 
an even parity function of y and square integrable (therefore normalizable) . Why is 
it then that this solution has been casted away from textbooks? The weaker answer 
would be because it does not provide quantized energy values, which is known to exist in 
any confined quantum system. The stronger answer is however that the function ipv(y) 
violates the boundary condition ^><(0) = ^>(0) for any non-integer u, which must be 
obeyed by the even parity wave functions ip(y) when g = 0. We are now about to see 
that ip' u (0 + ) = —ipl(0~). It is the latter property of ip u (y) which allows the solution of 
the quantum problem ^ with finite g. 

We have now gathered all the information needed to find the solutions of the 
eigenvalue problem ([8]). Since the Hamiltonian ^ is invariant over a parity the 
transformation, their eigenstates ip g (y) are either even or odd parity states. In the 
case of odd states we have ip g (0) = and therefore they do not see the presence of the 
delta-function at the origin. Thus, the odd parity wave functions ipg dd {y) are the states 
ip n {y) of the ordinary harmonic oscillator, with v = n = 1,3,5,... and the eigenvalues 
are e = n + | . To latter result we call the "trivial solution" of the eigenproblem (J8 ) . 

The solution of the even parity eigenfunctions ipg Ven {y) is not as simple, since these 
states feel the presence of the delta function at the origin. We need to find now the 
boundary condition the function w(y) must obey at the origin. Proceeding as we did at 
the introductory section, we integrate Eq. (Tsl) around y = obtaining: 



w>(0 + ) - w'JO-) = 2gw(0) . (15) 



Eq. (15) enable us to find the quantized energies of the even parity eigenstates we 
are seeking. Thus, the correct wave function for ^g Ven (y) is tp v (y) and not ip n (y) with 
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Figure 1. (color online) Graphical solution of Eq. (18 1, for several positive and 
negative values of g. The numerical values of v are given in Table [2] The eigenvalues 
for the energy E v are given by v + | with v the interceptions of the graphs with the 
x-axis (dashed lines refer to negative g, whereas solid lines refer to positive g). 



0, 2, 4, . . ., as in the case g = 0, for the latter wave function violates the boundary 



n 

condition (15). Using the results 0,/] 
lim .U(-v/2,l/2,y 2 ) = 

x— >0+ 

lim U'{-v/2,l/2,y 2 )-. 

x— >0+ 



7T 



T(l/2 - u/2) 



r(i - u/2) ' 

the eigenvalues associated with even parity eigenstates of Eq. 
numerical solution of the transcendent equation 

m-u/2) 



(16) 
(17) 

6 1) are given by the 



F(u) = u-g- 



0, 



T(l/2-u/2) 

which follows from the boundary condition (15). In Fig. Q we give the graphical 
solution of Eq. (18) for for g = ±0.25, g = ±1.0, g = ±2.5, and g = ±5.0, and in 
Table [2] the corresponding numerical values of u, for the first five even-eigenstates. As 
expected, the effect of the potential is to shift the eigenenergies of the even-states of the 

|| Here, the second identity can be derived from the first one using the fact that 9, g] U'(a, b, z) = 
-aU(a + 1,6+ l,z). 
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ordinary harmonic oscillator up or down in energy for positive and negative values of g, 
respectively. This effect is stronger for the low-lying eigenvalues (as we can anticipate 
from perturbation theory), and shifts the eigenenergies of the states if)„(x) toward those 
of their lower or higher neighboring odd-states, depending on the signal of g. This 
behavior is plotted in Figure |2j In the problem we are dealing with, and contrary to the 
simple case of the ordinary harmonic oscillator, if g < there is also a negative energy 
eigenvalue, as we could have anticipated from the solution of the attractive 5— function 
potential we have described in the introduction to this article. The absolute value of 
this negative energy state increases with the strength of the 5— function potential and, 
in the limit g — > — oo, the confinement imposed by the harmonic potential becomes 
irrelevant and the wave function transforms into the bound state given by Eq. [3] and 
with the same eigenenergy. Indeed, using Stirling's formula it is easy to prove 

that as z — > oo, T(z + \)/Y{z) — » y/z. For z = 1 + ||^| this implies that as g — > — oo, 
g ~ a/2|z/|. Since E ~ — the former result implies that 

rn 21 1 

E v = -g -hu = - — , 



which is the energy value we have obtain before for the simple case of an isolated 
attractive 6— function potential. 



g^O 


-0.25 


0.25 


-1.0 


1.0 


-2.5 


2.5 


-5.0 


5.0 


0.0 


-0.1557 


0.1281 


-0.8424 


0.3927 


-3.5865 


0.6434 


-12.9900 


0.7961 


2.0 


1.9288 


2.0693 


1.7208 


2.2546 


1.4285 


2.5042 


1.2305 


2.7003 


4.0 


3.9469 


4.0525 


3.7912 


4.2002 


3.5420 


4.4274 


3.3227 


4.6364 


6.0 


5.9558 


6.0439 


5.8258 


6.1699 


5.6051 


6.3772 


5.3833 


6.5887 


8.0 


7.9614 


8.0384 


7.8473 


8.1501 


7.6473 


8.3412 


7.4285 


8.5509 



Table 1. Eigenvalues v associated to the even parity eigenstates ipv(y) for several 
values of g; leftmost column corresponds to the case where g — > 0, which coincide 
with the harmonic oscillator even parity eigenvalues; the ensuing columns stand for 
g = ±0.25,g = ±1.0,3 = ±2.5 and g = ±5.0. In between each pair of eigenvalues we 
also have the states associated with the odd parity eigenstates, which have eigenvalues 
v = 1, 3, 5, 7, ... . The associated eigenenergies are given by e v = v + I 

Finally, in Figure [3] we plot in solid lines \ip n (y)\ 2 as function of y for n = 1,2,3. 
In the panels (a), (b), and (c) of the same figure we also plot, in dashed lines, \ip n (y)\ 2 
for the wave functions that would correspond to the harmonic oscillator with n = 2 for 
different positive and negative values of g. As g increases in positive (negative) values, 
the corresponding value of v approaches v = 3 [y = 2) making the absolute square of 
the wave function, |?/v(x)| 2 , looking like the absolute square of the wave function of its 

Note that in here, the harmonic oscillator frequency appears just as a by pass between the a and g 
coupling constants: at the end, the harmonic potential, too swallow compared with the delta, plays no 
role in obtain result. 
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:>0 — e<0 
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Figure 2. (color online) Graphical solution of Eq. (18 1, for several positive and 
negarive values of g. The numerical values of v are given in Table [2] 



odd-parity state neighbor, \ip n ±i(x)\ 2 . This doe not mean however that the two types 
of wave functions are the same, since they refer to orthogonal eigenstates. To make this 
point evident, we plot both type of states in panel (d) of Fig. [3j 

It is clear from Fig. [3] the behaviour ijj' u (0 + ) ^ (the number of nodes defines 
the order of the state). When g ^> 1, pairs of states (odd and even) of the harmonic 
oscillator with a delta-function become quasi-degenerated. Indeed, in this regime the 
dip of the wave function ip u (y) at y = approaches zero, but looking at Fig. [3] it is 
seen by the naked eye that the two functions are orthogonal (one is even and the other 
is odd; this is not self-evident from the density probability graphs). The enhancement 
of the curvature of the wave function around y = leads to an increase of the kinetic 
energy and therefore to an increase of energy of the even-parity eigenstates. 



3. Conclusions 

We have discussed the solution of the Schrodinger equation for the one-dimensional 
harmonic potential with a Dirac delta-function at the origin. The odd-parity eigenstates 
are given by the wave functions of the ordinary harmonic oscillator. This is obvious, 
since the states of the latter system are zero at the origin and therefore do not feel the 
presence of the delta-function. For the even parity states the solution is non-trivial. 
We have shown the existence of a solution of the differential equation of the harmonic 
oscillator that does not blow up at infinity as e x2 ^ 2 for non-integer values of v. As is well 
known, this solution is never mentioned in quantum mechanics textbooks for a good 
reason, but unfortunately that reason is, as far as we know, never discussed. Here we 
have shown that reason lies in the fact that its derivative at x = is finite, violating 
the boundary conditions imposed in the simple harmonic oscillator problem, that is, 



without the 5— function at the origin. Nevertheless, the wave function (14) is the one 



we need for solving the problem (pj) . In our work we have computed the eigenvalues and 
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(a) |g| = 1.0 (b) |g 2.5 




(c) |g| = 10.0 (d) fel =5.0 




Figure 3. (color online) Panels (a), (b) and (c): absolute square values of the wave 
functions of the harmonic oscillator (solid lines) for (from top to bottom) v = n = 3, 2, 1 
(plotted, respectively, in blue, black and red) and for tp v (x) (dashed lines) for v 
corresponding to values of g that varies from 1 to 10, for positive and negative 
values (top and bottom curves respectively). These plots show that the ^(x)! 2 
curves approach the neighbors odd states. However, as it is shown in panel (d), the 
wave functions are quite different since the ip v {x) are symmetrical with respect of x, 
presenting always a kink at x = 0. 



eigenfunctions of the even parity states and made, at the same time, a little excursion to 
the zoo of special functions using only the famous Handbook of Mathematical Functions 

Note added: After the submission of this article, the work by Busch et al. [5] was 
brought to our attention; the latter work elaborates on top of another paper by R. 
K. Janev and Z. Marie [TU] and, although focused on the three dimensional harmonic 
oscillator with a delta-function at the origin, one of the figures in Ref. [8] (Fig. 2) has 
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the same information as our Fig. [2j albeit presented in different form. Further more, 
the energy eigenvalues of the zero angular momentum states are given by an equation 
identical to our Eq. (18) but with 1/2 replaced by 3/2 (for obvious reasons). 
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